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PARTIAL DERIVATIVES OF MATRICES REPRESENTING RIGID BODY ROTATIONS 
By Samuel Pines, AMA, and B. F. Cockrell 


SUMMARY AND INTRODUCTION 


I*. 
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wT 

% 




The method presented for computing the partial derivatives of 
matrices representing rigid body rotations is applied here to the 
formulation of an RTCC processor presently under development for the 
lunar landing mission. This processor will determine the lunar module 
attitude with respect to the local vertical coordinate system by process" ng 
rendezvous radar shaft and trunnion angles with a weighted least squares 
filter. As will be seen, the necessary partials in final form are rather 
simple, but derivation would be cumbersome without the method derived 
here. 


In fact, to derive the required partials for the above mentioned 
processor with the standard method would require inversion of three 
3_by-3 matrices, the partial derivative of a 3-by-3 matrix with respect 
to its rotation angle, the multiplications of six 3-by-3 matrices, and 
a vector dot product . 

With the method given here, the derivation requires the definition 
of two unit vectors, a multiplication of one 3-by-l» matrix by a vector, 
and a dot product . 

The procedure described in this note can have wide application in 
all areas of dynamics involving rigid body rotations both inside and 
outside the realm of space technology. Examples include solving for 
tracking station location errors by processing angles, determining 
the final position of a rigid body after undergoing a finite rotation, 
determining inertial platform alignment errors with angle data, and 
determining any angular state which is a function of some angular 
observation. 
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RIGID ROTATION MATRICES IN THREE SPACE 


The differential equation for a rotation of a vector R, wich 

A 

respect to the variable o , about a unit vector, N, is given in any 
text on physics as 

A 

A * 

~ = N x R (1) 


The integration of this equation for a finite o will represent the rigid 

A A 

rotation of the vector R about the fixed unit vector N through the angle a. 


The vector cross product, in three space, is equivalent to a matrix 
transformation. If we take 


then 



"l 

*> 

N = 

n 2 

. |n| - i 


_ n 3_ 



m — 

r i 


R = 

r 2 



_ r 3 

> 


N x R 


n 2 r^ - n^ r 2 


n 3 r l " n l r 3 


n l r 2 


n 


2 



This may also be written as 


0 -n^ n 2 


"r/1 



1 

n_ 0 -n. 


r _ 

3 1 


2 

-n„ n., 0 


r _ 

L 2 1 


L 3_ 


( 2 ) 


N x £ 


(3) 
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If we write 


N * R = A(N)R 


then we can define the notation 


N * = A(N) 

where the matrix A is given by 


A ( N ) = n. 


(Ua) 


Using this notation, the differential equation for the rotating 
vector may be written as 

A 

'Nn * A 

I^AWR (5) 

dO 

This is a linear differential equation, and a solution may be 
obtained in terms of the unit characteristic solution 

AAA 

R = C(N, o) R ( 6 ) 


where the matrix C(N, o) satisfies the following conditions: 

8C( ” » - g ^ - A(N) C(N, o) (7) 

with the initial conditions given by the unit matrix in three space, i.e.. 


C(N, 0) = I. 


(Ta) 


The solution of equation (7) is given by 

„,: T \ oA(N) T 

C(N, o) = e I. 


Expanding the exponential matrix function, we have 

2 

A A fc- A A 


k 

a .k, 


C(N, a) = I, + oA(N) + 77 A (N) + . . . ♦ 77 A (N) 


( 9 ) 
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It car* be shown that the matrix A(N) satisfies the following 
characteristic matrix equation 

A A 

A 3 (N) = -A(N) 

A 

The solution of C(N, a) may be written as 

C(N, o) = 1 ^ + (l - ros o)A 2 (N) ♦ sin o A(N) 
This may also be written as 

% * A A 

C(N, o) = I^ ♦ (l - cos o) (N x) (N x) + sin a N * 

mJ 

A 

The inverse of C(N, o) is given by changing a to -o . 

C" 1 (N, o) * C(N, -o) 

It can be easily shown that this matrix is invertable by the 
following discussion: 


From equation (8), 


C(N, a) 


a A(N) 

e 


I 3 + o A(N) 


,[ q a(n) ] 2 

2 ! + 


+ [g A(N) n ] 
n! 


If we write A(N) in the classical diagonal canonical form 

T" 1 A(N)T = A 


or 

A (N ) = T A T" 1 

then [ A ( N ) ] 2 1 (T A T" 1 ) (TAT -1 ) = T A 2 T" 1 

In general [ A ( N ) ] n = T A n T” 1 

We now rewrite equation (8a) 


C(N, o) 


o T A T' 1 
e 


I + o(T A T" 1 ) 


o T A 2 T 
2 ! 



o T 

+ — 


Pre-and postmultiplying by T _1 and T, respectively 

T -1 C(N, o)T = I + oA+ — ^7- + . . . + 

d • 

or 


oA 


n 


n! 


T” 1 C(N, a )T = e 


o A 


( 10 ) 

( 11 ) 

(11a) 

(8a) 


(8b) 
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Now this is the classical canonical form of the matrix C ( N , a). 

A G A 

Therefore, the eigenvalues of C(N, o) are e , where A are the eigenvalues 

A A 

of A(N). It is easily seen that the eigenvalues of A(n) are 0, ± i. 

[see equation ( 10 ) ] 

The determinant of C(N, o) is then (e°) (e ia ) (e -10 ) = 1 and the 
matrix is invertable. 


We now prove that C(N, o) is an orthogonal transformation. We have 



T T 
C A 


( 12 ) 


T T 

But A is skew symmetric and A = -A where the superscript T denotes 
the transpose. It follows that 


We also have 


|£- = -C T A 
3o 


C" 1 C = I 


The derivative of equation (lU) yields 


fc.c-.c-c-.* 


= -C- 1 AC 


or 


£- C- 1 
3o C 


= C” 1 A 


It follows that 


.“1 


which proves that C is orthogonal. 

An interesting and useful consequence of equation (7) is 


SC( ^> a) C' 1 (N, o) = A(N) = H 


3o 


( 13 ) 

(lU) 

(15) 

(16) 

( 17 ) 


X 


( 18 ) 
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Moreover, if we have a product of two transformations, such as, 

E = D(M, S) C(N, o) 

the partial derivative of E with respect to the generic variable y 
given by 


(19) 

is 


|i . c(Ni 0) + D(M> fi) j 


( 20 ) 


Postmultiplying with the transpose of E, we have 


fy E = If B(M) + TT D(M ’ B) A(N) D? (M » 

A A 

where A(M) = M < 

We now prove that 

D(M, B) A(N) D T (M, B) = A(L) = L * 

where 

A A A 

L = D(M, B) N 


T 

The triple matrix product DAD is given by 


d ll 

d 12 

d 13 


0 

" n 3 

n 2 


’ d u 

d 21 

d 3l" 

d 21 

d 22 

d 23 


n 3 

0 

- n l 


d 12 

d 22 

d 32 

_ d 31 

d 32 

d 33 < _ 


_- n 2 

n l 

0 


_ d 13 

d 23 

d 33_ 


(21) 

(21a) 


( 22 ) 


(22a) 


(23) 


Carrying out the multiplication fct the first row of B we have 
b ll = V d 12 d 13 ' d 12 d 13 ) * n 2 (d ll d l3 ’ d n 


- n, (d. 


+ n 3 (d ll d 12 * d ll d l 2 > 


'12 - “1 '“13 d 22 - d 13 d 2 3 > + n 2< d u d^ - d„ d,J 


21 13 


+ n 3^ d 12 d 21 " d ll d 22^ 




b 13 = n l (d 13 d 32 " d 12 d 33 } + n 2 (d ll d 33 “ d i3 d 31 ) 
+ n 3 (d 12 d,^ - d ix d ?2 ) 


( 21 *) 
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Since the matrix D(M, 8) is orthogonal we have that 


*\ 


d 31 

Z 

d 12 

d 23 

- d 13 

d 22 

d 32 

S 

d 21 

d 13 

- d n 

d 23 

d 33 

s 

d n 

d 22 

- d 12 

d 21 

d 21 

s 

d 13 

d 32 

- d 12 

d 33 

d 22 

z 

d ll 

d 33 

d 13 

d 31 

d 23 

s 

d 12 

d 31 

- d n 

d 32. 


( 25 ; 


Consequently , 





b 12 = ’ d 31 n l ~ d 32 n 2 ' d 33 n 3 = ' 3 ) 
b 13 = d 21 n l + d 22 n 2 ' d 23 n 3 “ *2 


(26) 


Similarly we may show 



It follows that 



J 


D(M, 8) A(N) D T (M, 8) 


L « 


(27) 


('.3) 


where 


D(M, 8) N * L 


(28a) 


The final result for the product transformation is riven by 


3E 



f-M 




(29) 


APPLICATION 


The utility of this formulation may be seer in the discussion below 
which derives the partial derivatives to be used in an Apollo RTCC 
processor. 

Just after the lunar landing, the orientation of the LM is recorded 
with respect to a mean fixed coordinate system. If this stored alignment 
changes (due to LM settling) and realignment cannot be made with the 
alignment optical telescope; a separate method for LM altitude 
determination must be available. The reference presents a method in 
which the LM body attitude with respect to the local vertical coordinate 
system is determined by processing rendezvous radar (RR) shaft and trunnion 
angles with a weighted least squares filter. This method requires the 
partial derivatives of the observations with respect to a set of state 
elements. This state will be defined as three positive rotations about 
the LM local vertical system axes . The three rotations are ordered as 
follows : 


1 . 

About 

the 

local vertical, 

V 

2. 

About 

the 

displaced 

east , 

V 

3. 

About 

the 

displaced 

north 

’ “3 


The transformation from local vertical to LM body coordinates is 
then defined: 


cos 

a 3 

sin 

o' 


’cos 

0 

-sin a 2 


"l 


0 0 

-sin 

*3 

cos 

0 


0 

1 

0 


0 

cos 

sin 

_ 0 


0 

1 . 


.sin a 2 

0 

cos a 2 _ 


.0 

-sin 

cos a^_ 


(30) 

* A 

where R g is a vector in the body coordinates and R is the same vector 
in the local vertical system. 





From the body vector the observations may be written as 

sin(T) = -Y 


tan(S) = ~ 

Z B 

where S and T are RR shaft and trunnion angles. 


(30) 




* 

& 


to* 

!> 

I 




Given the above definition, the least squares method requires the 
following partia-L derivatives 


30 

3 cT 


where 0g is S or T and i = 1, 2, 3. 


B 


This may be written as 


90 

9a 


B 


90 


B 


9R, 


(31) 


3R 

9a 


B 


(32) 


Wk 

1 


I 

v 


The partial*: of 0g with respect to R p follow from the definition in 
equation ( 31 ) - They are 


9S 2. 

"t - * b~rf) [cos s> °* - sIn s] 

dR_ 


B 


— = [-sin T sin S, -cos T, -sin T cos S] 
3R B 


(33) 


(3U) 


Now by application of equation (5), 

9R . . 

— - A( V r b 

A ^ A 

Here N. is a unit vector about which Rg is rotated through 


(35) 


o. (i = 1, 2, 3) 
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Inspection of the vector columns in the matrices of equation (30) 
define the N^'s. They are 

N = [-cos a,, cos a.j, sin cos a^, - sin a,>] 

N 0 * [-sin -cos a^, 0] 

N 3 ■ [0, 0, -1] 

From the geometry of the LM body axes, Rg may be written as: 

cos T sin S 
-sin T 
cos T cos S. 



(36) 


(37) 


using the definition (equation Ua) for A(N.) and applying equations 
(36) and (37) to equation (35) 


3R 


B 


-H 1 = A(N, ) R p 

3a^ 1 B 


3H . . 

— * A( V r b 


3R 


B 


~ = A (N- ) R 

3a 3 B 


"-sin 

a 2 

sin 

T ♦ sin cos 

-sin 

“2 

cos 

T sin S + cos a 

_-sin 

°3 

cos 

cos T sin S 

"-cos 

°3 

cos 

T cos S 

sin 

°3 

cos 

T cos S 

_ cos 

a 3 

cos 

T sin S ♦ sin a 

f-sin 

T 


1 


■cos T sin S 
0 


(38) 


(39) 


(Uo) 
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The dot products of these vectors (equations 38, 39, and Uo) with 
trie vectors (equations 33 and 3*0 as defined by equation (32) forms 
the required six partial derivatives. 


as as 



3R 


— B 

- — = sin a_ cos a_ -sin a_ cos S tan T -cos a 0 cos a, sin S 

3 2 2 2 3 

(ia) 


aq as 3R b 

= — — • — — * -cos a. -sin S sin a_ tan T 

3o 2 3R_ 3 °2 3 3 

D 


(U2) 


as as 


aR 


B 


aa 3 3 R b aa 3 


-tan T cos S 


(U3) 


aR 

aT aT b 

t — = — t— • ? — = sin a_ sin S -cos cos a_ cos S 

3 “l 3R fi 3a l 2 2 3 


(1*U) 


3T 

aT 

0 

3R b 

3c 2 * 

A 

3R B 

3o 2 

3T 

3T 

A 

3 °3 " 

3r b 

aa 3 


-sin cos S 


sin S 


(1*5) 


( U6) 


tan T 
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